Plasma-Induced Frequency Chirp of Intense Femtosecond Lasers 
and Its Role in Shaping High-Order Harmonic Spectral Lines 
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We investigate the self-phase modulation of intense femtosecond laser pulses propagating in an 
ionizing gas and its effects on collective properties of high-order harmonics generated in the medium. 
Plasmas produced in the medium are shown to induce a positive frequency chirp on the leading 
edge of the propagating laser pulse, which subsequently drives high harmonics to become positively 
chirped. In certain parameter regimes, the plasma-induced positive chirp can help to generate 
sharply peaked high harmonics, by compensating for the dynamically-induced negative chirp that 
is caused by the steep intensity profile of intense short laser pulses. 
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Recent remarkable progress in high-power femtosec- 
ond laser technology has provided a novel opportunity to 
investigate high-order harmonic generation processes in 
an unprecedented high-intensity, ultrashort pulse regime. 
While necessary for generating high-order harmonics 
with wavelength extremely short compared to the laser 
wavelength , the use of lasers in this regime has some 
other important merits: the harmonic conversion effi- 
ciency is dramatically enhanced I^Q, and the harmonic 
pulse duration can be reduced to a subfemtosecond time 
scale U|. 

As the laser pulse evolves into the high-intensity, ul- 
trashort regime, high harmonic emission begins to ex- 
hibit complicated spectral features, such as harmonic line 
broadening, blueshifting, splitting, and smearing [pHlO[. 
All of these features are absent in the high harmonic spec- 
tra observed for lasers with weak intensity well below 
the ionization saturation intensity and long pulse dura- 
tion (longer than a few hundred femtoseconds) that show 
well-defined narrow harmonic peaks almost exactly at 
odd multiples of the fundamental laser frequency. The in- 
volved spectral features showing up for intense short laser 
pulses can be satisfactorily accounted for, at least at the 
single-atom level, in terms of decomposition of high har- 
monics into quantum-path components, and different be- 
haviors of the quantum-path components, in laser fields 
with rapidly- varying instantaneous intensity pO|-p^ . 

Since high-intensity laser pulses can ionize atoms 
producing plasmas that can affect propagating waves, 
a proper understanding of macroscopic high harmonic 
structure necessitates a detailed knowledge of the 
plasma-related propagation effect as well as the single- 
atom effect. It is known that an increase in electron den- 
sity during the laser pulse duration causes the refractive 
index of a medium to decrease with time, which leads to 
frequency upshifting and spectral broadening of the laser 
pulse p^-p^. Because the change in the driving laser 
spectrum should be reflected in the harmonic conversion. 



one might expect that harmonic spectral lines would like- 
wise be blueshifted and broadened. In reality, however, 
the harmonic spectrum exhibits behavior more compli- 
cated than this expectation, as demonstrated in this pa- 
per. As a matter of fact, there is another source that can 
affect the harmonic spectral line shape: the dynamically- 
induced harmonic chirp (dynamic chirp) that is brought 
about by a steep pulse envelope [ pl]jl^ . The observed 
spectral behavior of high harmonics can only be under- 
stood through inspecting the temporal variation of har- 
monic frequency caused by the plasma-induced change in 
the laser frequency, in close connection with the dynamic 
chirp. 

In this paper, we elucidate how the plasma effect mod- 
ulates the spectrum of an intense femtosecond laser pulse, 
and discuss its subsequent influence and the effect of laser 
focusing on macroscopic high harmonic spectra. It is 
shown that the plasmas induce a positive frequency chirp 
on the leading edge of the laser pulse up to the point at 
which a maximum blueshift is attained, and a negative 
chirp on the remaining part of the laser pulse. Depending 
on the relative amount of the plasma-induced chirp com- 
pared to the dynamic chirp, not only broadening but also 
narrowing can occur in high harmonic spectral lines. We 
demonstrate these using a one-dimensional (ID) model, 
in which the atomic response to the laser is calculated 
from the ID Schrodinger equation, and propagations of 
the laser and harmonic fields are considered in ID space 
along the propagation axis. 

In order to see how the laser field Ei is affected by 
the plasmas (whose effects are dominated by electrons) 
produced in the medium, we begin by finding a solution 
of the ID wave equation: 
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where ujp{x,t) = u)o[Ne{x,t)/NcrY^'^ is the local plasma 
frequency, and uiq is the laser frequency. The critical 
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plasma density Ncr is given in Gaussian units by Ncr = 
TOeajQ/47re^, where rUe is the electron mass. To calcu- 
late the electron density Ne{x, t), we use the ADK model 
pof , and consider sequential tunneling ionization up to as 
high stages of ionization as needed, neglecting collisional 
ionization that is of little significance in the parameter 
regions in which high harmonic generation experiments 
are commonly carried out. At the present gas pressures 
(< 100 Torr) much lower than 1 atm, the energy loss and 
temporal broadening of the laser are negligible jl^; thus, 
we may ignore the amplitude modulation of the laser 
field. Assuming the Gaussian incident pulse Ei{x,t) = 
Eoexp[-{2\n2/At'^){t - x/c)"^ - iwoit - x/c)] that is a 
solution to Eq. (1) in free space with At being the full 
width at half maximum (FWHM) of the pulse, we may 
then write the solution of Eq. (1) in the medium as 

Ei{x,t) = Soexp [-{2\ii2/At'^){t - x/cf ~ iujot 

+1 J\{x',t-\x-x'\/c)^dx'], (2) 

where n{x,t) = [1 — Nf,{x,t)/Ncr]^^'^ is the refractive in- 
dex of the medium. This expression turns out to be a 
good approximate solution of Eq. (1) under the condi- 
tions Ne/Ncr < 1 and \c~'^ J'' {dn/dt)dx'\ < 1, which 
are satisfied in the parameter regions considered in this 
paper. 

To confirm that Eq. (2) indeed closely approximates 
the exact solution, we present in Fig. 1 some typical spec- 
tra of 30- fs laser pulses at the exit of the medium (Ne gas) 
calculated from Eq. (2), along with those obtained from 
direct numerical calculations of Eq. (1). Use of the ex- 
plicit expression in Eq. (2) of course significantly reduces 
the computational time. Upon comparison, it is obvi- 
ous that the approximate solutions presented in Fig. 1 (a) 
agree well with the exact numerical results in Fig. 1(b) 
in the parameter regions considered. It can be seen that, 
as the gas density increases and/or as the laser intensity 
increases, the spectrum shifts toward a higher-frequency 
region and becomes broader. 

More detailed features of the laser pulse passing 
through the medium can be revealed with the help of 
the Wigner distribution function , that allows a view 
of temporal variation of the laser spectrum. The Wigner 
distributions, calculated under the same conditions as in 
Fig. 1, are displayed in Fig. 2. It can be observed that, 
owing to the phase modulation induced by plasmas pro- 
duced in the medium, the laser frequency increases with 
time (becomes positively chirped) in the leading edge, 
and then decreases (becomes negatively chirped) back 
to the original frequency in the remaining part of the 
pulse. We note that at the moment when the plasma- 
induced chirp changes sign from positive to negative, 
the production rate of electrons reaches its maximum. 
This moment, at which the laser experiences a maximum 
blueshift, comes earlier in time as the laser intensity in- 



creases, as can be seen from Fig. 2. 

When focused laser beams are used, as in usual high 
harmonic generation experiments, due regard should also 
be paid to the focusing effect that can change the am- 
plitude and phase of the laser field. In fact, in three- 
dimensional (3D) simulations this focusing effect could 
have been automatically considered. Apparently, how- 
ever, the ID wave equation in Eq. (1) cannot deal with 
this effect. Nevertheless, via Eq. (2) the focusing effect 
can be taken into account along the propagation axis as 
follows: 



Eif{x,t) = Ei{x,t)f{x)exp[i<^f{x)], 



(3) 



where f{x) = [1 + 4(a;/6)2]-i/2 and <i>f{x) = 
— tan~^(2a;/6) represent, respectively, the amplitude and 
phase changes due to the focusing , and b is the con- 
focal parameter. The laser field given in the form of Eq. 
(3) now suitably describes the phase modulation induced 
by plasmas and the focusing effect along the axis, and 
can be used for discussing the phase matching issue of 
high harmonics generated on the propagation axis. 

We next discuss the propagation of high harmonic 
fields. Since the change in the refractive index in the 
presence of plasmas is much smaller for higher-frequency 
waves, we may ignore here the plasma effect. Then 
the Green function method enables us to write a solu- 
tion of the 3D wave equation V^Eh — c""^ {d"^ Eh / dt^) = 
4TTC~^{d^P/dt^) in the following integral form: 
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where t' is the retarded time defined by i' = t — |r — r'|/c, 
and the laser-induced polarization P is given in terms of 
the gas density A'o and the atomic dipole moment d by 
P(r', t') = Nod{r', t'). Equation (4) clearly indicates that 
the harmonic field Eh at t is determined by the coherent 
sum of the dipole accelerations of atoms in the medium 
calculated at the retarded time t' . 

In our ID calculations, a medium lies on the propaga- 
tion axis, and the integration in Eq. (4) is performed in 
practice over the ID space along the axis. The medium, 
a Ne gas of length I = 700 /im at ~ 28 Torr, is uniformly 
discretized by 200 points, and at each point the dipole 
acceleration is calculated by numerically solving the ID 
Schrodinger equation for an atom in the laser field Ei / 
given in Eq. (3), with At = 30 fs, A = 800 nm, 6 = 4 mm, 
and / = 1 X 10^^ W/cm^ at the entrance of the medium. 
After suitably weighted by a constant factor to yield the 
correct gas density, the results are then added according 
to Eq. (4) to give Eh- 

The high harmonic spectrum calculated in the above 
way is presented in Fig. 3, which shows the respective 
roles played by the laser focusing and plasmas in the 
formation of macroscopic harmonic fields. In Fig. 3(b), 
we neglect the plasma effect by setting n{x,t) = 1 to 
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concentrate only on the focusing effect. Both the focus- 
ing and plasma effects are fully considered in Fig. 3(c), 
and for comparison we present in Fig. 3(a) a single-atom 
spectrum calculated for an atom located at the entrance 
of the medium. Whereas the single-atom spectrum in 
Fig. 3(a) is smeared with a complicated structure in the 
plateau region, the macroscopic harmonic spectrum in 
Fig. 3(b) exhibits discrete harmonic peaks. This is be- 
cause the intensity-sensitive harmonic components (the 
long quantum-path component and multiple-recoUision 
components), which give rise to the complicated struc- 
ture, are suppressed due to their poor phase-matching 
conditions, leaving only the short quantum-path compo- 
nent jl^JI^. Here we emphasize that the variation in 
the laser intensity caused by the focusing as represented 
by f^ix) is mainly responsible for this cleaning up of 
the harmonic spectrum. At very low gas densities, we 
may observe a harmonic spectrum similar to that shown 
in Fig. 3(b). However, at the gas density used in Fig. 
3(c), a further change in the harmonic spectrum is caused 
by the plasma effect: harmonics of low order (43rd and 
below) undergo spectral broadening, while higher-order 
harmonics get sharpened. 

The plasma-induced harmonic line broadening and 
narrowing are detailed in Figs. 4 and 5, respectively. 
The mechanism that shapes high harmonic lines can 
be explained by making a comparison between the 
plasma-induced chirp and the dynamically-induced chirp 
P, pl]p^ . Considering that only the leading edge of the 
laser pulse, where the depletion of neutral atoms is not se- 
vere and the electron density is low, is important for the 
phase-matched harmonic generation, we focus only on 
this part of the laser pulse. The leading edge of the laser 
pulse, in which the laser frequency increases with time 
due to the plasma effect (Fig. 2), tends to generate posi- 
tively chirped high harmonics. This effect is opposed by 
the dynamic single-atom effect that, in the present case, 
forces high harmonics to become negatively chirped . 
Depending on system parameters and harmonic order, 
one effect matches or dominates over the other. For in- 
stance, for the low-order harmonics in Fig. 4, the dynamic 
negative chirp is very small [Fig. 4(b)], and the plasma 
effect dominantly affects the time-frequency characteris- 
tics of high harmonics in such a way as to make high 
harmonics become positively chirped [Fig. 4(d)], leading 
to spectral line broadening [Fig. 4(c)]. In contrast, the 
dynamic negative chirp of higher-order harmonics in Fig. 
5(b) is substantially large, and the plasma-induced posi- 
tive chirp more or less cancels out the dynamic chirp [Fig. 
5(d)], resulting in harmonic line narrowing [Fig. 5(c)]. In 
both cases of Figs. 4 and 5, the plasma effect blueshifts 
the central frequencies of high harmonics. 

In conclusion, wc have investigated the plasma-induced 
phase modulation of intense femtosecond laser pulses and 
its effects on macroscopic high harmonic spectra, using a 
simple but accurate ID model. It has been shown that 



the plasmas produced in the medium induce a positive 
frequency chirp on the leading edge of the laser pulses. 
Strikingly, the plasma-induced chirp, which broadens the 
laser spectrum, can lead not only to broadening but also 
to narrowing of harmonic spectral lines. The underlying 
mechanism has been explained by comparing the plasma- 
induced positive chirp with the dynamically-induced neg- 
ative chirp. In stark contrast to the widely reported find- 
ing (in regimes well below saturation) that high harmon- 
ics are negatively chirped |2lj2^j2^, this study clearly 
demonstrates that, under laser fields of high intensity 
and gas pressures not too low, high harmonics can be- 
come free of chirp, or even become positively chirped, 
thanks to the plasma effect. This suggests a way to con- 
trol the chirp and spectral line shape of high harmonics 
using the plasma effect. 
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FIG. 1. (a) Laser spectra obtained from Eq. (2) for a 
30-fs (FWHM) laser pulse of wavelength A = 800 nm after 
propagating through a Ne gas medium of length I — 700 ^im. 
The laser intensity and gas density are given by J = 1 x 10 
W/cm^ and A^o = 1 x 10 



1 X 10'^ W/cm^ and A^o = 3 x 10' 



cm"^ (~ 28 Torr) (thin solid line); 

" cm"^ (~ 85 Torr) 
(dotted line); J = 3 x 10^^ W/cm^ and Ao = 1 x 10^* cm'^' 
(dashed line). The spectrum of the incident laser pulse is 
drawn by a thick solid line, (b) Laser spectra obtained by di- 
rect numerical calculations of Eq. (1) for the same parameters 
as in (a). 
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FIG. 2. Wigner distributions of the laser pulses specified 
in Fig. 1. In (a)-(c) [(d)-(f)], the same parameters are used as 
for the thin solid, dotted and dashed lines in Fig. 1(a) [Fig. 
1(b)], respectively. Only positive contour lines are shown. 
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FIG. 3. (a) Single-atom harmonic spectrum for a Ne atom 
placed at the entrance of the medium, (b) Macroscopic har- 
monic spectrum (in the far field, on the axis) considering only 
the effect of laser focusing, (c) Macroscopic harmonic spec- 
trum taking into account both the focusing and plasma ef- 
fects. A Ne gas medium of density A^o = 1 x 10^* cm~^ and 
length / = 700 fj,m is irradiated by a 30-fs (FWHM) Gaussian 
laser pulse of wavelength A = 800 nm. The medium is cen- 
tered at a; = 2 mm behind the focus, and because of focusing 
{b — 4 mm) the laser intensity decreases from / = 1 x 10^^ 
W/cm^ at the entrance to 1 = 0.7 x lO^'' W/cm^ at the exit 
of the medium. 
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FIG. 4. (a) Enlarged view (on a linear scale) of Fig. 3(b) 
between 20a;o and SOcJo (only the focusing effect is consid- 
ered), (b) Wigner distribution of the harmonics chosen in 

(a), (c) [(d)] Same as (a) [(b)] except that harmonics arc 
chosen from Fig. 3(c) (both the focusing and plasma effects 
are considered). Positive and negative values of the Wigner 
distribution are colored black and white, respectively. 
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FIG. 5. Same as Fig. 4 except that harmonics between 
TOwo and SOwo are chosen here. 



